Multipole Analysis of Radio-Frequency Reactions in Ultra-Cold Atoms 
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Using the multipole expansion we analyze photo induced reactions in an ultra-cold atomic gas 
composed of identical neutral bosons. While the Frank-Condon factor dominates the photo induced 
spin-flip reactions, we have found that for frozen-spin process where the atomic spins are conserved 
the reaction rate is governed by the monopole r 2 and the quadrupole terms. Consequently the 
dependence of the frozen-spin reaction rate on the photon wave number k acquires an extra k 4 factor 
in comparison to the spin-flip process. Comparing the relative strength of the r 2 and quadrupole 
modes in dimer photoassociation we predict that the mutual importance of these two modes changes 
with temperature and scattering length. 
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Radio-Frequency (rf) association of molecules in ultra- 
cold atomic traps is an invaluable experimental tool for 
studying the universal properties of few-body systems 
[3, In such experiments pj-flfl] rf induced dimer or 
trimer formation leads to enhanced atom loss rate due 
to larger probability for three-body recombination. The 
measured atom loss rate as a function of the rf field fre- 
quency incorporates the data about the various molecular 
thresholds and structures, and can be used to calibrate 
the two-body scattering length and to measure the three- 
body binding energy. 

Photo-integration of molecules, where resonant rf radi- 
ation is used to stimulate photon emission, is closely re- 
lated, through time reversal, to electro-weak reactions in 
light nuclei. In fact, rf photorcactions inducing an atomic 
spin-flip are governed by the Frank-Condon factor. This 
factor is the atomic analog of the Fermi operator, which 
at low momentum transfer describes the charged vector- 
current neutrino induced isospin-flip reactions in nuclear 
systems. The nuclear analog to frozen-spin rf reactions is 
photoabsorption, where an incoming photon dissociates 
the nucleus. For example, photoreactions within the a- 
clustcr model of nuclei such as 12 C, 16 are governed by 
the same physical mechanisms studied here. Therefore, 
studying rf photoreactions in ultra-cold atomic systems 
can be instrumental for understanding the nuclear dy- 
namics, the complicated nuclear continuum, and for pre- 
dicting unmcasurable electro-weak cross-sections of as- 
trophysical importance. 

Previous analysis of these rf experiments [Ill - [l4j , have 
relied on the Franck-Condon factor estimating the transi- 
tion rates. While this approach is appropriate for describ- 
ing spin-flip reactions, it is not suitable for frozen-spin 
reactions. In [15[ we have studied the quadrupole con- 
tribution to the frozen-spin response of a bound bosonic 
trimer. In this article we present a complete multipole 
analysis of the rf association process forming a bound 
molecule composed of N identical bosons. We show that 
the spin-flip and frozen-spin processes differ by their op- 
erator structure and by the de-excitation modes that con- 
tribute to the photoassociation rate. We apply our re- 
sults to study frozen-spin dimer formation. We focus on 



bosonic alkali-metal atoms, in which a single s-shell va- 
lence electron is coupled to an half-integer spin nucleus 
to form a boson. This is the case for 7 Li, 23 Na, 39 K, 41 K, 
85 Rb, 87 Rb, and 133 Cs. Extension of our results to other 
atoms and fermionic isotopes is rather simple but lays 
outside the confines of the current manuscript. 

In ultra-cold atomic traps, the trapped particles are 
subject to an external magnetic field B = Bz applied 
to control the inter-atomic interaction, i.e. the scatter- 
ing length, through a Feshbach resonance. Consequently, 
the energy of the system is dominated by the interactions 
of the valence electron spin S — 1 /2 with the external and 
the nuclear magnetic fields. The Zeeman interaction with 
the external magnetic field reads —gfisS ■ B, where q is 
the electron g- factor, and /is is the Bohr magneton jla ]. 
The hypcrfinc spin-nucleus contact interaction fnohfl ■ S 
couples the electron with the nuclear spin I where u>hf 
is the hyperfine splitting frequency. When both interac- 
tions are comparable, the total spin F = I + S is not a 
good quantum number, however its projection nip, is. 

To be specific, we consider / = 3 /2, the case for 7 Li, 
23 Na, 39 K, 41 K, and 87 Rb. The atomic spin state com- 
position in the \ms,nii) basis, 

\m F ) = sin6» mF |i/2,m F - i/ 2 ) + cos 6 mF |-i/2, m F + 1/2) 
is given by 
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(1) 

for the low energy branch states. Here 77 = g^sB /2huj} l f 
is the ratio between the Zeeman and the hyper-fine split- 
ting. 

In photoassociation experiments, an rf radiation of few 
MHz is applied to the ultra-cold atomic gas. When the 
photon energy matches, stimulated emission occurs and 
a molecule is formed. Fermi's golden rule dictates the 
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transition rate of such process to be, 

r ^f = T EE k/'MWI 2 ^ ~ E f- fa*) ■ ( 2 ) 

where averages on the appropriate initial states and 
sums on the final states. Specifically, \i) is the con- 
tinuum state, kp is the emitted photon momentum and 
polarization, and |/) is the formed A^-body bound state. 

The interaction Hamiltonian Hj — — - J dxJ(x) ■ 
A{x) between an atomic system and an electromag- 
netic (EM) radiation field A is defined by the current 
operator J composed of convection and magnetization 
terms, J(x) = J c {x) + cV x fi(x). For a system 
of neutral atoms interacting with long wavelength ra- 
diation J c « 0, and the interaction between the ra- 
diation field and the system occurs through the mag- 
netization density e^i(x) = gnB^2iSi5(x — r*). Us- 
ing box normalization of volume f2, the EM field reads 

M*) = £ fc , P \fM;^ P (*U tk - x + h.c) where p = 1, 2 
are the photon linear polarizations, Uk is its frequency 
and k is its momentum. 

For an atomic system we can assume that the ini- 
tial and final matter wave functions can be written 
as a product of spin and configuration space terms, 
* = ®lm l Xm f , where $ LA / L = $lm l (ri, r 2 , . . . ta) 
is a symmetric iV-particle wave function with angu- 
lar momentum quantum numbers LAIl and xm f = 
S [|771f(1))|771f(2)) . . . Im^tyl))] is the symmetrized spin 
wave function with magnetic quantum number Mp = 
^2jinp(j). Using this factorization, The transition ma- 
trix element in ^ takes the form 



calculated to give, 
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The spin operator can be written in spherical nota- 
tion, S ■ (fc x e kp ) = £ A (-) A 5-A ■ ( k x £kp)\ where 
A = 0, ±1. To calculate this quantity the laboratory z 
axis is defined by the static magnetic field. The rotation 
from the laboratory reference frame (x, y, z), to the pho- 
ton reference frame, (e, fc x e, fc) is represented by Eulcr 
angles (a, /?, 7). Then, (fc x e) ■ z = — sin 7 sin /3, and 
|(fc x e)±i| 2 = (1 -sin 2 7sin 2 /3)/2. 

For rf radiation the photon wavelength is much larger 
then the typical dimension of the system R, therefore 
kR <C 1 and the lowest order in kR dominates the in- 
teraction. The first contribution depends on the photon 
energy versus the energy difference between adjacent Zee- 
man states. 

Using Eq. (TT]), the spin matrix elements can be easily 



{sin 6 m > p cos 9 mp A = +1 
-±cos26> mF A = 
sin 6 mF cos 6 m ' p A = -1 
(4) 

If the photon can induce a Zeeman state change, i.e. spin- 
flip, the leading contribution comes at order k. Energy is 
delivered to the system through the spin matrix element 
and we can approximate e lx ~ 1, to get 
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The factor A^ 2 results from utilizing the symmetry of the 
wave function and Sj,\ is the magnetic moment operator 
of any specific particle j . The last term on the rhs of Eq. 
([5]) is just the Frank-Condon factor. 

In frozen spin reactions A = and the rf photon does 
not affect the spin structure of the system. Therefore 

M' F = Mp, \xm f ) = \x\i F ) an d the transition matrix 
element can be written as 



(f,kp\Hi\i) = - igp, B 
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where (Sq) = (Xm f \Si,o\Xm f ) is the average single par- 
ticle magnetic moment. 

In the long wavelength limit, the leading contributions 
to the photoreaction process are 
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where Yi m arc the spherical harmonics. The difference 
between the four operators in is clear. The first op- 
erator is proportional to 1 and stands for elastic interac- 
tion. For photon emission this kind of interaction defies 
energy conservation. The second term is just the dipolc 
operator, which for identical particles is proportional to 
the center of mass and hence does not affect relative mo- 
tion of the atoms. Therefore the only contributions to 
such excitation arises from the r 2 , and the quadrupolc 
terms. Consequently the transition probability scales as 
fc 5 . Summing over the initial and final magnetic numbers 
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Ml, M' l , the transition matrix element reads 



M L M' r 



2il 



(8) 



Here, like in ([5]), we have used the symmetry of the wave 
function and rj is the coordinate of any of the parti- 
cles j = 1, ... N. Comparing equations flSJ) and ([5]) the 
difference between spin-flip and frozen-spin reactions is 
evident. To begin with, the frozen-spin acquires an ex- 
tra fc 4 factor suppressing the reaction rate. Then, we see 
that whereas the leading configuration space spin-flip op- 
erator is proportional to the unit operator, like the Fermi 
operator in nuclear physics, for frozen-spin reactions we 
identify two competing modes of order k 2 , one is the r 2 
operator and the other is just the quadrupole operator. 

Once established the operator structure of the transi- 
tion rate we would like to focus our attention on frozen- 
spin photoassociation reactions trying to estimate the rel- 
ative importance of the two reaction modes. To this end 
we consider the case of dimer formation, which under 
certain simplifying assumptions can be calculated ana- 
lytically. 

In the universal regime, where the scattering length 
is much larger than the potential range, one can ap- 
proximate the wave functions by their asymptotic be- 
havior. The dimer bound state with binding energy Eb 
is an s-wave state, and <£>Q (r) = Yoo Vs( r )/ r where 
r = T2 — v\ is the relative coordinate. <^s(^) — V2~Ke~ Kr 
and k = \J m\EB\/ti 2 - For the initial (continuum) states, 
the asymptotic scattering solution is used, i.e. & LM (r) = 
Y LM (f>L(r)/r, where <f> L (r) = y/2/ Rqr[j L {qr) cos 8 L - 
yL(qr) sin 5 l], 8l is the phase shift, ji and i/l are the 
spherical Bessel functions of the first and second kind, 
q is the relative momentum, and the wave function is 
normalized in a sphere of radius R. 

Given the L' = dimer ground state, the r 2 term 
corresponds to L = s-wave formation, while the 
quadrupole term to L = 2 rf-wave formation. Using 
|(0||r ||0)| 2 = 1/4* and |(2||F 2 ||0)| 2 = 5/4n, 
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= J Q <p* B (r)4>o(r)r 2 dr, and same for Id replacing 
2. The transformation to the relative coordinate 
r m ((yj) adds an extra factor of 1/16 in (jHJ. 

For the s-wave mode, the radial integration yields 



Here I. 

(j> by i 



(q 2 + re 2 ) 3 



(3k 2 — q 2 ) cos So (3q 2 — n 2 ) sin 8i 

(10) 



At low scattering energy, the s-wave phase shift takes the 
form qcot^o = — 1 /a B + r e ffq 2 /2 where a s is the scatter- 
ing length and r e // is the effective range. The dimer's 
binding energy is connected to these parameters by the 



relation [li 



y as +r eff n 2 /2. 



For the c?-wave mode, the radial integration yields, 



cos 82 




(11) 



For short-range potentials, 8l 
^5 



q2L+i ag g — y 0, thus 
^2 ~ {a-ciq) 5 where ad is the d-wave "scattering length". 
We note that for the Van der Waals potential the correct 
threshold behavior is, 8 2 ~ ±(a d g) 4 [17j . As long as ad 
is of the order of the effective range it is safe to assume 
8 2 « 0. 

The relative contribution of the s, d modes to the dimer 
formation is displayed in Fig. [U where the last term in 
parenthesis on the rhs of Eq. (J5J is presented normalized, 
along with the s and d components. From the figure it 
can be seen that the s-wave association is peaked around 
q = k/2, while the d-wave association is peaked at q = k. 
For small a s /r e ff ratio the d-wave is the dominant mode 
but as we approach the unitary limit a s /r e ff — > 00 the 
s-wave becomes dominant. We note however, that as the 
a s /r e ff ratio become smaller finite range corrections, ig- 
nored here, become more important. In photoassociation 
experiments angular analysis is impossible and therefore 
the different energy dependence of the two dimer forma- 
tion modes, Fig. [I] can be used to distinguish between 
them. As it stands, the large overlap between the two 
modes will certainly complicate such an attempt. 

The average of the initial states and the sum over 
the final states J^f are the last components needed for 
evaluating the dimer association rate, Eq. The ini- 
tial two-body state \i) = \qLMi) describes two atoms 
in the continuum with relative momentum q and angular 
momentum L, Ml. The average on these states takes the 
form 



E = E E 



L=0 M L = -L 



dqP(qLM L ) , 



(12) 



where P{qLMi,) is the probability of finding an atomic 
pair in an internal state \qLMi) due to thermal distri- 
bution. Wc assume that the system is in thermal equi- 
librium with temperature T > Tbec higher then the 
condensation temperature. The partition function of a 
pair confined in a sphere of radius R is given by 
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dq J2(2L + 1)(R 
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dq 



4 





1.0 r 




0.8 - 




0.6 - 




U.4 - 






itio 


n 9 
u.z - 


jj 
s 

CO 


o.o L 


p-< 




+^ 




TJ 


0.8 - 






s 


0.6 - 


o 


0.4 - 








0.2 - 




0.0 L 





FIG. 1. (Color online) The normalized 2-body transition ma- 
trix element, Eq. ((8]), as a function of the relative momen- 
tum q/n. The sum (black), the r 2 (red, peaked at q = k/2) 
and quadrupole (blue, peaked at q = k) terms are given for 
a / r eff — 2 (upper panel) and a/r e // = 200 (lower panel). 



Given q, finite R imposes a cutoff L max w qR on L. For 
L = 0, 2 and large enough i? we can safely assume that 
P(qLM L ) = P(q) = £S e -/s»V/ro 

The sum ^ j contains all possible dimer quantum num- 
bers and all possible photons weighted by the emission 
function, 



E 



(2tt) 3 



]T fdk(l + N kp ), 

t i 7i r/ *' 



(14) 
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where N kp is the number of photons with momentum k 
and polarization p in the initial state. The stimulating rf 
radiation is a narrow distribution centered at some k r f. 
Therefore N kp w Nk rfPrf S(k — k r f)5 PtPr} /Q, where fc r / 
is the momentum of the stimulating rf field and is its 
polarization. 

Substituting Eqs. ([8]), ([12]), ([14]) into Fermi's golden 
rule ([3]), the dimer formation rate is given by 



2 H 2 B ck^. f mN Kf 



Pr-f 



16(2ir) 2 H 2 q 
\(So)\ 2 



r 7 sin /3( ^7, + 152^ 



(15) 



The relative momentum g is connected to the photon 
wave number through energy conservation h 2 q 2 /m = 
Eb + Hck r f. 

The relative importance of s, d modes shifts with tem- 
perature. This point is demonstrated in Figs. [5] and [3] 
In Fig. [5]we present the s, d rates to the photoassociation 
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FIG. 2. (Color online) Normalized dimer photoassociation 
rate as a function of rf photon frequency. The total (black), 
the r 2 (red, low energy peak) and the quadrupole (blue) rates 
are presented for 7 Li with a s = 1000ao at T = 5/iK (lower 
panel) and at T = 25/iK (upper panel) . 



of 7 Li dimers at a — 1000ao, ao is the Bohr radius, for 
T = 5/iK and T = 25/iK. The relative contribution of the 
r and quadrupole to the resonance peak is displayed in 
Fig. [3] It is evident that at low temperatures the s mode 
is dominant. The imprtance of the d mode grows with 
the ratio UbT/Eb- Therefore we can conclude that for 
small ksT/Es values the photoassociation is an s-wave 
process while for large values it is a d-wave process. 



Summing up, in this work we have studied the photo 
response of ultra-cold atomic gases. Comparing spin- 
flip and frozen-spin experiments we have found that the 
reaction mechanism is quite different in these two pro- 
cesses. While the spin-flip reaction is dominated by the 
well known Frank-Condon factor, for frozen-spin process 
the reaction rate is governed by the two competing r 2 
and quadrupole modes. Comparing the relative strength 
of these modes in dimer photoassociation we predict that 
the mutual importance of these two modes changes with 
temperature and scattering length. The implications of 
these results on photoassociation experiments and on 
trimer formation rates is yet to be studied. 
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